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Abstract 

We study the differential branching ratio, branching ratio and the CP-violating asymmetry for 
the exclusive Bd —>■ (??, ry')£^£" decays in the standard model. We deduce the Bd —>■ [rj, rj') form 
factors from the form factors of 5 ^ tt available in the literature, by using the SU{'i) f symmetry. 
We observe that these decay modes, which are within the reach of forthcoming B -factories, are 
very promising to observe CP-violation. 
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1 Introduction 



The decays of B-meson are very promising for investigating the Standard Model (SM) and search- 
ing for the new physics beyond it. Among these B-decays, the rare semileptonic ones have at- 
tracted much attention for a long time, since they offer the most direct methods to determine 
the weak mixing angles and Cabibbo-Kobayashi-Maskawa (CKM) matrix elements. These de- 
cays can also be very useful to test the various new physics scenarios like the two Higgs doublet 
models (2HDM), minimal supersymmetric standard model (MSSM)| 1 1. etc. 

From the experimental side, there is an impressive effort to search for B-decays, in B-factories 
such as Belle, BaBar, LHC-B. For example, CLEO Collaboration reports for the branching ratios 
(BRs) of the 7r-£+iy and S° p^i+u decays as 

BR{B^^TT~£+u) = (1.8±0.4±0.3±0.2) X 10""^ 
BR{B^^p-tu) = (2.57 ± 0.29l|]j^ ± 0.41) X 10"^ (1) 

From these results, the value of the CKM matrix element \ Vuh\ = 3.25±0.14^q 29 ±0.55 has been 
determined |2|. Recently, the BR of the inclusive B Xsi^i^ decay has been also reported by 
Belle Collaboration [3 1 ; 

BR{B Xsl^r) = (6.1 ± lAt\\) X 10^^ , (2) 

which is very close to the value predicted by the SM 0. 

In this paper, we investigate the B^ — > rj^'^ decay modes within the SM. It is well 
known that the inclusive rare decays are more difficult to measure, although they are theoretically 
cleaner than the exclusive ones. This motivates the study of exclusive decays, but their theoretical 
investigation requires the additional knowledge of decay form factors, i.e. the matrix elements of 
the effective Hamiltonian between initial B and the final meson states. The nonperturbative sector 
of QCD is used in order to determine these form factors. Two of the form factors, and /_, 
necessary for B^ ^ -q £+^^decay have been calculated very recently, in the framework of light- 
cone QCD sum rules [SJ. However, we do not have a precise calculation on the remaining form 
factor, fx for B^ — > rj^'^ decay yet. Therefore, in this work, we choose to deduce the form 
factors of Bd rj^'^ transition from the form factors of B —> n using the SU (3) f symmetry. The 
form factors of B ^ it have been calculated in the light-cone constituent quark model (LCQM) 
|6l0 and also in the QCD sum rules method (QCDSR) |8|; and in this paper, we will give our 
numerical results using both of these approaches. Let us mension that the B —>■ K hadronic 
matrix elements computed in LCQM and QCDSR have been used to evaluate the semileptonic 
rate of the B —>■ Kt^i~ decay mode [9, 10|. Compared to the recently measured value of 
BR{B K£+£-) = {0.75toil ± 0.09) x IQ-^ by Belle Collaboration Ol and also BaBar 
Collaboration 1 12], we see that QCDSR predicts a better result. 

In this work, we also calculate the CP asymmetry in the B^ rj^'^ decay, which is 
induced by the b — > ci^'''^" transition at the quark level. For b — > si'^i^ transition, the matrix 
element contains the terms that receive contributions from ti, cc and uu loops, which are propor- 
tional to the combination of = VtbV^*, = VcbVcs and = KibKTs' respectively. Smallness 
of in comparison with and ^t, together with the unitarity of the CKM matrix elements, bring 
about the consequence that matrix element for the b —>■ si'^i' decay involves only one indepen- 
dent CKM factor ^t, so that the CP violation in this channel is suppressed in the SM 1131 IT4l . 
However, for b d£~^£~ decay, all the CKM factors r]t = VfbV^*^, rjc = VcbVcd ^"^^ ^« = ^ubyud 
are at the same order in the SM so that they can induce a CP violating asymmetry between the 
decay rates of the reactions b —>■ di^£~ and b —>■ d£^£~ |15|. So, b d£^£^ decay seems 
to be suitable for establishing CP violation in B mesons. On the other hand, it should be noted 
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that the detection of the b —>■ d£^£^ decay will probably be more difficult in the presence of a 
much stronger decay b si^£^ and this would make the corresponding exclusive decay chan- 
nels more preferable in search of CP violation. In this context, the exclusive —>■ (vr, p) £'^£~ , 
and —>■ £'^£~ decays have been extensively studied in the SM I16IIT71 and beyond |[T8l - l22l . 

The paper is organized as follows: In section |2l first the effective Hamiltonian is presented 
and the form factors are defined. Then, the basic formulas of the differential branching ratio 
dBR/ds, branching ratio BR and the CP violating asymmetry Acp for Bd tj^'^ decays 
are introduced. Section|3lis devoted to the numerical analysis and discussion. 

2 Effective Hamiltonian and Form Factors 



The leading order QCD corrected effective Hamiltonian, which is induced by the corresponding 
quark level process b —>■ d£~^£~, is given by li23J - l.26J : 



AC ( ^° 



i(^)[on/^)-Oi(/^)] 



(3) 



where 



Au 



VubV. 



ud 



VtbV, 



(4) 



td 



using the unitarity of the CKM matrix i.e. V^bV^^ + Kb^J^ = explicit forms of 

the operators Oj can be found in refs. l23l l24l . In Eq.©, Cj(/i) are the Wilson coefficients 
calculated at a renormalization point n and their evolution from the higher scale /x = my/ down 
to the low-energy scale /i = rrih is described by the renormalization group equation. For Cj^^{fi) 
this calculation is performed in refs. ll27ll28l upto next to leading order. The value of Cio{mi,) to 
the leading logarithmic approximation can be found e.g. in |.23„2_6 ,| . The terms that are the source 
of the CP violation are given by the following, which have a perturbative part and a part coming 
from long distance (LD) effects due to conversion of the real cc into lepton pair £~^£~ : 



C'g"i^i) = C^'\^i)+Yresonis) , 



(5) 



where 



Cg + h{u, s)[3Ci{fi) + C2ifi) + 3C3(m) + C^ifi) + 3C75(^) + ^(m) 



+ A„(3Ci + C72 
1 



-h{l, s) (ACsifi) + 4C4(m) + 3C5(m) + Ceifi)) 



-h{0,s) [C^ifi) + 3C4(/i) + A„(6Ci(/i) + 2C2(m))] 
y 



(6) 



and 



Y 



.is) 



a 



3 ^ 



7rT{Vi^£+£-)mv, 



em v;=^, m%s- mv, + imyTv, 



X [(3Ci(/i) + C2(/i) + 3C3(^) + C4(/x) + 3C5(/i) + C^{^x)) 

+ \u{3Cl{^,) + C2{^i))]. 



(V) 
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In Eq.@, s = /m?B where q is the momentum transfer, u = ^ and the functions h{u, s) arise 
from one loop contributions of the four-quark operators Oi — Oq and are given by 

w N 8 , mfo 8 , 8 4 

h(u,s) = — m inu-\ — -H — v (8) 

^ ' ^ 9 9 27 9^ 



2 -, ,„ hi 



ivr , fory = ^<l 



2 arctan / . foryEE^>l, 

7 N 8 8 , m{, 4 , 4 . 

27 9 /i 9 9 

The phenomenological parameter k in Eq. Q is taken as 2.3 (see e.g., IT51l ). 

Neglecting the mass of the d quark, the effective short distance Hamiltonian for the b 
dt^i~ decay in Eq.© leads to the QCD corrected matrix element: 

M = ^I^VM\clff{mb) d7^(l - 75)6^r^ + C7io(mb) d^^{l - 75)bhV 

(10) 

Next we proceed to calculate the BR?, of the Bd rj^'^ decays. The necessary ma- 
trix elements to do this are < r?W(p^(,))|(i7^(l - 75)&|-B(pb) >, < r]'^'\p^(,))\di(Tf,^q^{l + 
75)6|i?(pB) > and < '?/^'HPr;('))l'^(l + l5)b\B{pB) >■ The first two of these matrix elements can 
be written in terms of the form factors in the following way 

<V^'\Pr,io)\d-/i,{l-75)b\B{pB)> = /+(g')(ps+P^w)/. + r(5')g/., (11) 

< V^'Hrr^oMa^.^q^'il + -f5)b\B{pB) > = [(pb + P^w)/.'?^ " q^.{ml - m^^o^)]f^{q^Xi2) 

where ps and p^{,) denote the four momentum vectors of B and r/^'^ -mesons, respectively. fv{q'^) 
is sometimes written as fv{q^) = fr/im-B + ^"^o )• 

To find < ri^'\p^(^,))\d{l + j^)b\B{pB) >, we multiply both sides of Eq. dTTl with q^ and 
then use the equation of motion. Neglecting the mass of the d-quark, we get 

< 7?W(p^(,))|d(l + 75)b\B{pB) >= ^[/+(g2)(m| - m^o^) + f- {q^)q% (13) 

As pointed out in sec.l, although the form factors /+ and /_ for i? — > ry decay have been 
calculated in the framework of the light-cone QCD sum rules in fSl, we do not have a precise 
calculation of the other form factor fy in the literature yet. However, the form factors of B^. rj^'^ 
transition can be related to those of 5 ^ vr through the SU{3)f symmetry l29ll30l . In addition, 
the authors of emphasize that, their results coincide with the ones that are calculated using the 
SU{3)f symmetry. Therefore, we choose to deduce the form factors necessary in this work from 
the B TT transition using the SU{3)f symmetry. For rj — r]' mixing, we adopt the following 
scheme EIlEa, 

\t] > = cos (j)\r]q > - sin (plijs >, 

\r]' > = sin (p\r]q > + cos (p\r]s >, (14) 
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where \r]q >= {uu + dd)l\f2, \t]^ >= ss, and (/) = 39.3 is the fitted mixing angle 1311 . Hence, 
the relation between the form factors are written as follows: 



F^d^^^q^) = cos 0F^^^(g2), 

^BW(g2^ = sin (15) 

For S — > vr, we use the results calculated in two different frameworks: In the LCQM, the form 
factors are parametrized in the following pole forms (6l0l 



J ^1 ) r ^2 X2.35 '/ W ; . .2.30 ' 

Uq') = - , _"7x2.31 • (16) 

[} ~ Km) 

However in the QCDSR approach, they are given by |8| 

, 2n 0.305 ^ , 2. 0.305 

1 _ 1.294- + 0.206(4-)2 1 - 0.266^- - 0.752(4-)2 

= , ^ (17, 



0.05 



1-1.284- + 0.193(^)2^ ' 



from which / can be calculated through the relation: 

f- = {ml-mlMo-f+)/q\ (18) 

Using the above matrix elements, we find the amplitudes governing the Bd rj^'^ decays 
as follows: 



(19) 



where 



A = C'/ff+ -2mBC'/^f., 
G = Cio/+, 

D = Cio(/+ + r). (20) 

Using Eq.(fT9t and performing summation over final lepton polarization, we get for the double 
differential decay rates: 

= ^IW/-|VA.|m|A(l-zV)|AP 

+ (m|A(l - z^v'^) + 16 r ml) |G|^ + A s mj \D\'^ 

+ 4.m}{l-r - s) Re[GD*]\ , (21) 



Here s = /rn\, r = m^(,)/m^, v = \/l — t = mj/m'^^, A = r^ + (s — 1)^ — 2r(s + l), 
and z = cos 6, where 6 is the angle between the three-momentum of the lepton and that of 
the B -meson in the center of mass frame of the dileptons i^£^ . After integrating over the angle 
variable we find 

^^B-.r,iO ^2 2 



ds 2i07r5 



E^\Vti,V:fm%^vA, (22) 



where 



A = lm|A(3-.2)(|^|2 + |G|2) + M(12.. + A)|Gp 
S 6s 

+ Amjs\D\^ + Amj{l-r-s)Re[GD*]. (23) 

We now consider the CP violating asymmetry, Acp, between the rj^'^ £^£~a.nd 

jj(') £+£" decays, which is defined as follows: 

''''^ ' T{Bd ^ r?W £+£-) + T{Bd ^ r/W £+i~) ' ^ ' 

Using this definition we calculate the Acp as: 

_[ 

J{A-H{s))ds 

where 



H{s) = I /+m|(3 - v^)X Im A„ (im 6 C'j" fr — /+ {Im Ci ^2)] -(26) 

3 \^ niB + J 

In calculating this expression, we use the following parametrizations: 

Cf s f 1 + A„ {2 , (27) 

K- ^-'-'^-fr' +OiX^ (28) 

(1 - p)^ + T]^ 

3 Numerical Results and Discussion 

In this section we present the numerical results of our calculations related to B^ — > rj^'^ £'^£~ {£ = 
e, /X, r) decays, for four different sets of parameter choice of the form factors and the updated 
fits of the Wolf ens tein parameters fT^ . which are summarized in Tabled The total BRs are 
collected in Table |2l We have also evaluated the average values of CP asymmetry < Acp > in 
Bd V^'^ £"^^~decays for the above sets of parameters, and our results are displayed in Table 
|3l In both tables, the values in the paranthesis are the corresponding quantities calculated without 
including the long distance effects. We observe that the results of < Acp > is very sensitive to 
the choice of four different sets of parameters for r channel, while they are very close to each 
other for // channel. 

The input parameters and the initial values of the Wilson coefficients we used in our numerical 
analysis are as follows: 

ruB = 5.28 GeV , nib = 4.8 GeV , = 1.4 GeV , = 1.78 GeV, 

= 0.105 GeF, \VtbVt*a\ = 0.01, = 0.547 GeF, m^, = 0.958 GeF, 
Ci = -0.245, G2 = 1.107, C3 = 0.011, G4 = -0.026, C5 = 0.007, 
Ce = -0.0314, Gj^^ = -0.315, Gg = 4.220, Gw = -4.619. (29) 
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(p; v) 


Form factors 


set-1 


(0.3; 0.34) 


LCQM 


set-2 


(0.15; 0.34) 


LCQM 


set-3 


(0.3; 0.34) 


QCDSR 


set-4 


(0.15; 0.34) 


QCDSR 



Table 1: List of the values for the Wolfenstein parameters and the form factors of the transition 
B ^ TV calculated in the light-cone constituent quark model (LCQM) |l6l|7| and light-cone QCD sum 
rule approach (QCDSR) [i8j|. 



10« ■ BR 


e 


setl 


set2 


set3 


set4 




T 


0.331 


0.313 


0.687 


0.659 






(0.324) 


(0.314) 


(0.695) 


(0.677) 


V 




2.704 


2.511 


3.704 


3.468 






(2.119) 


(2.063) 


(3.049) 


(2.966) 




e 


2.713 


2.520 


3.716 


3.479 






(2.127) 


(2.371) 


(3.059) 


(2.976) 




r 


0.092 


0.087 


0.153 


0.146 






(0.086) 


(0.083) 


(0.147) 


(0.144) 


1]' 


/i 


1.363 


1.268 


1.779 


1.666 






(1.033) 


(1.010) 


(1.395) 


(1.365) 




e 


1.369 


1.273 


1.786 


1.674 






(1.038) 


(1.015) 


(1.402) 


(1.372) 



Table 2: The SM predictions for the integrated branching ratios for i = t, fx,e of the B^ rj^'Hi 
decay with (without) the long-distance effects. 

There are five possible resonances in the cc system that can contribute to the decay under 
consideration and to calculate their contributions, we need to divide the integration region for s 
into three parts for i = e, fj. so that we have 4m|/m^ < s < {m^^ — 0.02)^/m^ and (m^^ + 
Omf/m% < s < (m^2 - Omf/ml and (m^^ + 0.02) Vm-| < s < (m^ - m^(,))V"^B> 
while for £ = T it takes the form given by 4m^/m^ < s < imip2 ~ 0.02)^/m^ and {m^^ + 
0.02)^/m^ < s < {tub — m^(/))^/m^ . Here, m^^ and m^,^ are the masses of the first and the 
second resonances, respectively. 

In Fig. O and Fig. we present the dependence of the BR on the invariant mass of 
dileptons, s, for the ^ rj r+r" and B^ ^ rj /i+/i^decays, respectively. We plot these 
graphs for the parameter set-1 and set-3 in Tabled represented by the dashed and the solid curves, 
respectively. The sharp peaks in the figures are due to the long distance contributions. As can be 
seen from these graphs, BR stands more for the parameter set-3. The same analysis above is 
made for B^ — > t]' t~^t~ and i?^ rj' decays in Fig. ^ and Fig. (|4]i, respectively 

Figs. ^ and ^ are devoted to the Acp{s) as a function of s for Bd —>■ 7?r+r^ and 
Bd ^ Tj P^A^ decays, respectively. In these figures, the small dashed (dotted dashed) and the 
solid (dashed) curves represent the Acp{s) for the parameter set-1 and set-3 with (without) long 
distance contributions. The dependence of Acp on s for the rj' channel is plotted in Fig. (Q and 
Fig. (Ell, for ^ = T and ^ = /U, respectively We see from these figures that for i = fi, Acp{s) is 
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10- < Anp > 




setl 


set2 


set3 


set4 




T 


1.291 


0.961 


2.271 


0.840 








I'D (S^yi 


(Tl SQ?"! 
I u.oy 1 I 


(Tl ^(SlTi 


V 


/i 


u.u^ ( 




u .uyz 


u.ozo 






fO 663) 


(0.484) 


(0.671) 


(0.490) 








u .^y u 


u.uyo 








fO 663) 

\ KJ • W \J '-.J J 


(0.484) 


(0.671) 


(0.490) 




r 


0.926 


0.693 


0.886 


0.656 








I'D '^1 CC\ 




yJAoo J 






0.578 


0.444 


0.593 


0.452 






(0.637) 


(0.464) 


(0.639) 


(0.465) 




e 


0.579 


0.444 


0.594 


0.452 






(0.638) 


(0.464) 


(0.640) 


(0.465) 



Table 3: The same as Table Q, but for < Acp >■ 



not very sensitive to the choice of the parameters set-1 or set-3, reaching up to 28% for the larger 
values of s for both the rj and -q' channels. However for £ = r case, Acp{s) gets slightly larger 
contribution from set-3 than set-1, but reaches at most 25% in the small-s region. We note that 
Acp{s) is positive for all values of s, except in some resonance regions. We also observe from 
table |3] that including the long-distance effects in calculating < Acp > changes the results only 
by 2 — 10% for ^ = ;U mode, but for £ = r, it becomes very sizable, 30 — 150%, depending on the 
sets of parameters used for (p; rj). 

In conclusion, we have analyzed the r]^'^ f+^^decays within the SM. We have found 

that, these decay modes have a significant Acp, especially for £ = t. Since calculated BRs 
of these decay modes are within the reach of forthcoming B-factories such as LHC-B, where 
approximately 6 x 10^^ Bd mesons are expected to be produced per year, we may hope that it can 
be measured in near future. 
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Figure 1: Differential branching ratio for B rjr^T^ decay as a function of s for the parameter 
set-1 and set-3, represented by the dashed and the solid curves, respectively. The sharp peaks in the 
figures are due to the long distance contributions. 
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Figure 3: The same as Fig.dll) but for the B ^ rj' t^t^ decay 
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s 

Figure 4: The same as Fig.([l]) but for the B ^ rj' ix'^fi~ decay 




s 



Figure 5: Acp{s) for B r]T^T~ decay for the parameter set-1 and set-3 with (without) long 
distance contributions, represented by the small dashed (dotted dashed) and the solid (dashed) curves, 
respectively. 



11 



1 r 




s 



Figure 6: The same as Fig.© but for the B ^ r] fi^fi decay 
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Figure 7: The same as Fig.© but for the B ^ rj' t^t^ decay 
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Figure 8: The same as Fig.© but for the B ^ r]' ix decay 
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